{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "Control valve design\n", "====================\n", "\n", "This is example 8.2 in Seborg, but worked a little differently to allow choice of $R$ and $C_{cv}$" ] }, { "cell_type": "code", "execution_count": 1, "metadata": { "collapsed": false }, "outputs": [], "source": [ "import numpy\n", "import scipy.optimize\n", "import matplotlib.pyplot as plt\n", "from ipywidgets import interact\n", "%matplotlib inline" ] }, { "cell_type": "code", "execution_count": 2, "metadata": { "collapsed": true }, "outputs": [], "source": [ "# Constant pump head\n", "DeltaPa = 40\n", "# Guess for q\n", "q0 = 100" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "The MEB reduces to quadratic form:\n", "\n", "$$ \\Delta P_a = \\Delta P_{hc} + \\Delta P_{v} $$\n", "$$ \\Delta P_a - a_{hc}q^2 - a_{v}q^2 = 0$$\n" ] }, { "cell_type": "code", "execution_count": 3, "metadata": { "collapsed": true }, "outputs": [], "source": [ "def MEBcoeffs(l, R, Ccv, characteristic='eqperc'):\n", " ahc = 30/200**2\n", " if characteristic == 'linear':\n", " fl = l\n", " elif characteristic == 'eqperc':\n", " fl = R**(l - 1)\n", " av = (1/(Ccv*fl))**2\n", " \n", " return [-ahc - av, 0, DeltaPa]" ] }, { "cell_type": "code", "execution_count": 4, "metadata": { "collapsed": true }, "outputs": [], "source": [ "def positive(x):\n", " return x[x>0][0]" ] }, { "cell_type": "code", "execution_count": 5, "metadata": { "collapsed": true }, "outputs": [], "source": [ "ls = numpy.linspace(0.01, 1)" ] }, { "cell_type": "code", "execution_count": 6, "metadata": { "collapsed": false }, "outputs": [], "source": [ "def curve(R, Ccv, characteristic):\n", " qs = [positive(numpy.roots(MEBcoeffs(l, R, Ccv, characteristic))) for l in ls]\n", " plt.plot(ls, qs)\n", " plt.plot([0, 1], [0, max(qs)])\n", " plt.xlabel('Valve lift (l)')\n", " plt.ylabel('q')\n", " plt.ylim([0, 250])" ] }, { "cell_type": "code", "execution_count": 7, "metadata": { "collapsed": false }, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYgAAAEKCAYAAAAIO8L1AAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMS4yLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvNQv5yAAAIABJREFUeJzt3Xl8VfWd//HXJ/tKIBACYV/CroJExH1BURG1zmirvxlXWmxHW22dTrX7bzp2nF+nnWnn12mro1U7jltXgmtFLW4oq2wBCYsQtgAJ2fd8549zEgIeIEDOvTfJ+/l43Me99+Tcez+HhLxzzvd7Psecc4iIiBwpLtoFiIhIbFJAiIhIIAWEiIgEUkCIiEggBYSIiARSQIiISKDQAsLMhpnZm2ZWZGbrzOxef/n3zWynma3yb3M6vOZBMys2s41mdkVYtYmIyPFZWOdBmNlgYLBzboWZZQLLgc8AnwWqnXP/esT6k4BngBlAHvA6MM451xJKgSIickyh7UE453Y751b4j6uAImDIMV5yHfCsc67BObcVKMYLCxERiYKESHyImY0EpgEfAOcB95jZrcAy4H7nXDleeCzp8LISAgLFzOYD8wHS09OnT5gwIdTaRUR6muXLl+93zuUcb73QA8LMMoDfAfc55yrN7BfADwDn3/8YuBOwgJd/6viXc+4R4BGAgoICt2zZsrBKFxHpkczsk86sF+osJjNLxAuHp51zvwdwzu11zrU451qBRzl0GKkEGNbh5UOBXWHWJyIiRxfmLCYDHgOKnHM/6bB8cIfVrgfW+o8XADeZWbKZjQLygQ/Dqk9ERI4tzENM5wG3AGvMbJW/7JvAzWY2Fe/w0TbgLgDn3Dozex5YDzQDd2sGk4hI9IQWEM65dwgeV3jpGK95CHgorJpERKTzdCa1iIgEUkCIiEggBYSIiARSQIiISCAFhIiIBFJAiIhIIAWEiIgEUkCIiEggBYSIiARSQIiISCAFhIiIBFJAiIhIIAWEiIgEUkCIiEggBYSIiARSQIiISCAFhIiIBFJAiIhIIAWEiIgEUkCIiEggBYSIiARSQIiISCAFhIiIBFJAiIhIIAWEiIgEUkCIiEggBYSIiARSQIiISCAFhIiIBFJAiIhIIAWEiIgEUkCIiEggBYSIiARSQIiISCAFhIiIBAotIMxsmJm9aWZFZrbOzO71l2eb2Z/NbJN/389fbmb2MzMrNrPVZnZmWLWJiMjxhbkH0Qzc75ybCMwE7jazScADwCLnXD6wyH8OcBWQ79/mA78IsTYRETmO0ALCObfbObfCf1wFFAFDgOuAJ/3VngQ+4z++DnjKeZYAfc1scFj1iYjIsUVkDMLMRgLTgA+AXOfcbvBCBBjorzYE2NHhZSX+siPfa76ZLTOzZfv27QuzbBGRXi30gDCzDOB3wH3OucpjrRqwzH1qgXOPOOcKnHMFOTk5XVWmiIgcIdSAMLNEvHB42jn3e3/x3rZDR/59qb+8BBjW4eVDgV1h1iciIkcX5iwmAx4DipxzP+nwpQXAbf7j24A/dVh+qz+baSZQ0XYoSkREIi8hxPc+D7gFWGNmq/xl3wQeBp43s3nAduBG/2svAXOAYqAWuCPE2kRE5DhCCwjn3DsEjysAzApY3wF3h1WPiIicGJ1JLSIigRQQIiISSAEhIiKBFBAiIhJIASEiIoEUECIiEkgBISIigRQQIiISSAEhIiKBFBAiIhJIASEiIoEUECIiEkgBISIigRQQIiISSAEhIiKBFBAiIhJIASEiIoEUECIiEkgBISIigRQQIiISSAEhIiKBFBAiIhJIASEiIoEUECIiEkgBISIigRQQIiISSAEhIiKBFBAiIhJIASEiIoEUECIiEkgBISIigRQQIiISSAEhIiKBFBAiIhIotIAws8fNrNTM1nZY9n0z22lmq/zbnA5fe9DMis1so5ldEVZdIiLSOWHuQTwBXBmw/N+cc1P920sAZjYJuAmY7L/mP80sPsTaRETkOEILCOfcYqCsk6tfBzzrnGtwzm0FioEZYdUmIiLHF40xiHvMbLV/CKqfv2wIsKPDOiX+sk8xs/lmtszMlu3bty/sWkVEeq1IB8QvgDHAVGA38GN/uQWs64LewDn3iHOuwDlXkJOTE06VIiIS2YBwzu11zrU451qBRzl0GKkEGNZh1aHArkjWJiIih4toQJjZ4A5PrwfaZjgtAG4ys2QzGwXkAx9GsjYRETlcQlhvbGbPABcDA8ysBPgecLGZTcU7fLQNuAvAObfOzJ4H1gPNwN3OuZawahMRkeMz5wIP9XcLBQUFbtmyZdEuQ0SkWzGz5c65guOtpzOpRUQkkAJCREQCKSBERCSQAkJERAIpIEREJJACQkREAikgREQkkAJCREQCKSBERCSQAkJERAIpIEREJJACQkREAikgREQkkAJCREQCKSBERCRQpy4YZGZfO9bXnXM/6ZpyREQkVnT2inIFwFl4lwYFuAZYDOwIoygREYm+zgbEAOBM51wVgJl9H3jBOff5sAoTEZHo6uwYxHCgscPzRmBkl1cjIiIxo7N7EL8BPjSzPwAOuB54MrSqREQk6joVEM65h8zsZeACf9EdzrmV4ZUlIiLR1tk9CJxzK4AVIdYiIiIxROdBiIhIIAWEiIgEUkCIiEggBYSIiARSQIiISCAFhIhId1K9j3ef+zG71vwl9I/q9DRXERGJkvJPYMNCXFEhbvsHnEcrSw58jrzTLgr1YxUQIiKxxjkoLYKiQthQCHvWALA3ZSzPNn8Gm3gN99z0mdDLUECIiMSC1lbYucwPhYVQtgUwGDaDlsv+kYe3juXRdTD/wtE8eNUEzCz0khQQIiLR0tIE2972Q+ElqN4DcYkw6kI498sw/moaUgfw5f9ZyWvr93L/5eO459KxEQkHUECIiERWYw0UL/L2Ej5+BeorIDENxl4GE6+B/NmQ2heA2sZm5j+xjHeK9/O9ayZxx3mjIlqqAkJEJGy1ZfDxq96ewuY3oLkOUvvBhLnebcwlkJh62Esqapu488mlrNxezo9uOJ0bC4ZFvGwFhIhIGCp3wYYXvVDY9g64FugzBM68xQuFEedBfPCv4D0V9dz2+Ids3V/Dz//PmVx12uAIF+8JLSDM7HFgLlDqnJviL8sGnsO72NA24LPOuXLzDqj9FJgD1AK3+91jRUS6j/3FULTAO3y0c7m3bMA4OO9emDgX8s6E44wfbNlXzS2PfcjB2kZ+fcdZnDd2QAQKDxbmHsQTwP8Hnuqw7AFgkXPuYTN7wH/+DeAqIN+/nQ38wr8XEYldzsHuVVC00AuFfRu85XnT4NLveGMKOeM7/XarSw5y+6+XYsCz88/htKFZ4dTdSaEFhHNusZmNPGLxdcDF/uMngbfwAuI64CnnnAOWmFlfMxvsnNsdVn0iIieltQW2v+/PPHoRKnaAxcOIc6HgTphwNWQNPeG3fbd4P/OfWkbftCR+M28Go3MyQij+xER6DCK37Ze+c263mQ30lw8BdnRYr8Rf9qmAMLP5wHyA4cOHh1utiAhAUz1secs7aW3jy1B7AOKTYcylcPEDMO4qSO9/0m//4urdfPW5VYwakM5T82aQ2yel62o/BbEySB10UM4FreicewR4BKCgoCBwHRGRU1ZfCZte8/YUil+HxmpI7gPjrvAGmcdeBsmn/lf+E+9u5f8uXM/04f147LazyEpL7ILiu0akA2Jv26EjMxsMlPrLS4COc7iGArsiXJuI9HbVpbDxJW9MYetfoKUR0gfCaTfAhGu8E9gSkrrko1pbHf/8chGPvr2Vyybm8h83TyM1Kb5L3rurRDogFgC3AQ/793/qsPweM3sWb3C6QuMPIhIR5dsODTJvXwI46DcSZsz3BpmHngVxXfuLu76phftf+IgXV+/m1nNG8L1rJhMfF5mzo09EmNNcn8EbkB5gZiXA9/CC4XkzmwdsB270V38Jb4prMd401zvCqktEejnnoHS9HwqHGuGROwUu+oY3HTV3ynGno56sg7WNzP/Ncj7cWsaDV01g/oWjI9Y640SFOYvp5qN8aVbAug64O6xaRKSXa2+Et8ALhvKteI3wzobZ/+SNKWSH38aipLyW23+9lO0HavnZzdO49oy80D/zVMTKILWISNdqbvQa4W1Y+OlGeOd9BcZfDZm5EStnTUkFdz65lIamFp6aN4OZo09+1lOkKCBEpOdoa4RXVOj1PmqogMR0yL/MG2QeNxtSIn/y2Str9/DV51aRnZ7E058/m3G5mRGv4WQoIESke6st87qiFi2EzYugud5rhDfx6I3wIsU5xy//soV/eWUDU4f15dFbC8jJTI5KLSdDASEi3U97I7wFsO3dDo3wbvOCYfi5R22EFymNza186w9reGF5CdeckcePbjidlMTYmsZ6PAoIEeke9m86dLW1k2yEFyllNY188b+9mUr3zsrnvsvyY3am0rEoIEQkNrU3wiv0Dh/t3+gtzzsTZn3XG1PIGRfdGgMUl1Yz78ml7K6o56c3TeW6qUOiXdJJU0CISOxoafYa4W1Y+OlGeGfNO+lGeJGyqGgv9z27iuTEOJ75wkymj+gX7ZJOiQJCRKKrrRFeUaHX5qKuDBJS/EZ4D8L4qyAtO9pVHlNrq+Pnbxbzk9c/ZnJeHx65pYC8vtEZGO9KCggRibz6Ctj059Ab4UVCTUMz9z//Ea+s28P104bwz391WrcbjD4aBYSIREZ7I7xC2PIXaG3yG+Hd6A0yj+y6RniR8smBGuY/tZxNpVV8++qJzDt/VLccjD4aBYSIhOdojfDOviu0RniRsvjjfXz5mZUAPHXn2ZyfH71Lg4ZFASEiXae9EZ4/82hvWyO807wL60yYC7mTY2Y66sloaXX8xxub+OmiTYwbmMmjtxYwvH9atMsKhQJCRE5NayuULPU6o0axEV4klNU0cu+zK3l7037+atoQ/un6KaQl9dxfoz13y0QkPG2N8NpmHlXv9Rrhjb7IO3Ft/JyINsKLhBXby7n76RUcqG7kh9efxs0zhvWo8YYgCggR6ZzGGm/GUdHCmGqEFzbnHE++t42HXioit08Kv/vSuZw2tOdtZxAFhIgcXXsjvELY/IbfCC/bG2CeOBdGXxy1RniRUFHXxDd/v4YX1+xm1oSB/OSzU2PqmtFhU0CIyOEqdnpnMW8o7NAIb2hMNcKLhOWflPGVZ1axp7Keb1w5gbsuHE1cDF4WNEw9/7ssIsfX1givqBB2rfCWDRgP59/nDTLnTevWM49OREur4xdvFfNvr28ir28KL3zxHM4c3r1bZpwsBYRIb+Qc7FrpnZ/QjRrhhW1PRT33PbeSJVvKuOaMPB66fgp9UnrPIaUjKSBEeou2RnhFhd4hpMoSrxHeyPPgrM/7jfC6b+fRU/Xn9Xv5+m8/orG5lR/dcDo3TB/a42cpHY8CQqQna6qHLW96ewlHNsK75JvdohFe2Krqm/inhUU8t2wHk/P68B83T2N0TvfoAxU2BYRIT9PeCG8BbHodmmogOctrhDfRb4SXlB7tKmPCki0H+PsXPmLXwTq+dPEY7rssn+SE7tn6IwwKCJGeoLrUn3m08FAjvIxcOP2z3bYRXpjqm1r411c38ti7WxmRncYLXzyH6SN6955UEAWESHfV1givqBB2fIDXCG8UzPyiN8g89CyIi4t2lTFndclBvvb8RxSXVnPLzBE8OGdCj26XcSr0ryLSXTgHe9cdmnnUAxvhham+qYWfLdrErxZvIScjmafunMGF43KiXVZMU0CIxLK2RnhFC7xgKN8GGAyfCbMf8mYe9ZBGeGH6YMsBHvz9Grbsr+HG6UP59tWTetUZ0SdLASESa5obYdviQzOPDmuEd58XChkDo11lt1BV38TDL2/g6Q+2Myw7lf+e1zOv2xAWBYRILGhvhFcIH7/WoRHe5V7fo/zLe2QjvDC9vn4v3/7jWkqr6vn8+aP42uxxGms4QfrXEomW2jLY+LJ36KgXNsILy86DdfygcD2vrNvD+NxMfnnLdKYO6xvtsrolBYRIJLU1witaAJ+8d6gR3vTbvUHm4ef0ikZ4YWhsbuWxd7bys0WbcDi+fsV4vnDBaJISNJPrZOknUSRs+z4+dLW1Xt4ILyzvbd7Pd/+0juLSamZPyuU7cycxLLtnXgY0khQQIl3tsEZ4hbD/Y2/5kOkw63veIaQB+dGtsYfYW1nPD18q4k+rdjEsO5XHby/g0gk960p20aSAEOkKLc2w/T1vL+FTjfC+0Osb4XW1usYWHn17C794azMtzvGVWfn83cVjSElUm4yupIAQOVntjfAKvcHm9kZ4s+DSb8G4K3t9I7yu1trqWPDRLv7llQ3srqhnzmmDeODKiQzvr8NJYYhKQJjZNqAKaAGanXMFZpYNPAeMBLYBn3XOlUejPpGjqq/wpqFuKFQjvAhb/kk5P1i4nlU7DjJlSB/+/XNTOXt0/2iX1aNFcw/iEufc/g7PHwAWOeceNrMH/OffiE5pIh20NcIrKoStiw81wjvjc94g88gL1AgvRFv2VfPjP3/Mi6t3MzAzmX+98Qz+atqQXnf5z2iIpUNM1wEX+4+fBN5CASHRUrb1UM8jNcKLit0Vdfxs0SaeX1ZCckIcX5mVz10XjiY9OZZ+bfVs0fqXdsBrZuaAXznnHgFynXO7AZxzu80ssJeAmc0H5gMMHz48UvVKT3dYI7xC2LvWWz7oNLj4Qe/w0cBJmo4aAeU1jfzyL5t54r1ttDrHLTNHcM+lYxmQkRzt0nqdaAXEec65XX4I/NnMNnT2hX6YPAJQUFDgwipQeoHWVij50L8EZ0AjvIlzod/IKBfZe1TWN/Hku9t4ZPEWqhubuX7aEL562TidzxBFUQkI59wu/77UzP4AzAD2mtlgf+9hMFAajdqkh2tvhFcIG16CmlK/Ed7FcP5XYfwcNcKLsIq6Jn797lYef2crlfXNXDYxl69fMZ7xgzKjXVqvF/GAMLN0IM45V+U/ng38I7AAuA142L//U6Rrkx6qodprhLdhoRrhxZCK2iYee3crv353K1X1zVw+KZd7Z+UzZYi+F7EiGnsQucAfzDuWmwD8j3PuFTNbCjxvZvOA7cCNUahNeoq2RnhFhd65Cs31kNYfJl3jDTKPvhgSU6JdZa+0v7qBJ97dxpPvbaOqoZkrJufylVn5TM5TMMSaiAeEc24LcEbA8gPArEjXIz1IRcmh6ahqhBdzPjlQw6Nvb+GFZSU0trRy1ZRBfPnSfCYO7hPt0uQo9L9Furf2RniFXv8jgJwJ3njCxLkweKpmHkXZ6pKD/OovW3h57W4S4uL46+lD+PwFoxmTkxHt0uQ4FBDSvTjndUQtWuiNKagRXkxqaXW8saGUx9/ZyvtbDpCZksBdF43hjnNHMrCPDu11FwoIiX3tjfAK/UZ4O/1GeOfDjPnezCM1wosJFXVNvLBsB0++v40dZXXkZaXwrTkTuWnGMDJTdA3o7kYBIbGpqQ42v+ntJXyqEd631QgvxmzaW8UT723j9yt2UtfUwoxR2XzzqolcPimXhHidcd5dKSAkdrQ1witaAMWLDjXCG3+lN8g8dpYa4cWQ+qYWXl23h6c/2M6HW8tISojjM1PzuO3ckZqR1EMoICS6qvbCxhe9MQU1wusWikurefbD7fxuRQnltU0Mz07jH64cz01nDSc7Xd+rnkQBIZHX3givEHZ8CDjIHg0zv+QNMg8pUCO8GFPT0MzLa/fw/LIdfLi1jIQ444rJg7h5xnDOHdNfnVV7KAWEhM85r/ld28wjNcLrFlpbHUu2HOC3K0p4Ze0eahtbGNk/jW9cOYEbpg8lJ1PN83o6BYSE46iN8M6BK37oXYJTjfBiUnFpNX9cuZM/rNzJzoN1ZCYncN3UPP76zKFMH9EPU5D3GgoI6TrNjd44woYOjfDik2DURWqEF+N2lNWycPVuCj/axfrdlcQZnJ+fwz9cOZ4rJg/StZ57KQWEnJq2RnhFhbDpNWiohKQMrwHehLmQPxtS1EohFu2trOelNV4orNh+EIBpw/vy3bmTuPr0weTqhLZeTwEhJ67mAHz8sjemsPkNaGnwG+Fdq0Z4Me6TAzW8um4Pr6zd0x4KEwf34R+uHM81p+fp2gtyGAWEdM5hjfDeBdcKWcOg4E5vkHnYTDXCi0HOOTbsqeK1dXt5Zd0einZXAjBlSB/+fvY4rpwyiLEDdd0FCab/0XJ0+zYeGmQ+rBHe19QIL4bVN7Xw/uYDLNqwlzeKStlVUY8ZFIzox7evnsgVkwdpT0E6RQEhh3RshFdUCAc2ecuHFMBl3/cOHw0YG80K5Sh2lNWyeNM+3txQyjvF+6lvaiUtKZ7zxw7g3svyuWTCQAZm6rCfnBgFRG/X0uwdMtqw8NON8M6+y5uO2icv2lXKEWoamnl/8wHe3rSPxZv2s3V/DQBD+6XyuYJhXDoxl7NHZWv2kZwSBURv1NYIr6jQG2yuK/ca4Y29DC79Doy7Qo3wYkxDcwurth/k/S0HeG/zAVZuL6epxZGaGM/M0dnces4ILsjPYUxOus5TkC6jgOgt6g5601CLCg81wkvJ8rqiqhFezGlsbmXNzgqWbDnA+5sPsOyTMuqbWjGDyXl9mHf+aC4cN4DpI/qRnKC9BAmHAqIna2+EVwhb3/Yb4Q2CM27yBplHXgDx6tEfC6obmlnxSTlLt5WxdFsZq3YcpL6pFYAJgzK5ecZwzhndn7NH9ScrTd8ziQwFRE9TtuVQzyM1wotJzjk+OVDLyh3lrNx+kBXby1m/q5JWB3EGk/OyuHnGcGaMzGbGqGz6Z6jnkUSHAqK769gIr6gQStd5ywedDpd80zt8NHCipqNGUVlNI2t2VrB6x0FW7jjIyu3llNc2AZCeFM8Zw/pyz6X5nDWyH9OG9yMjWf8tJTboJ7E7am3x9g7aWmYf/AQ1wosN5TWNrN1VwZqdFawp8e5LyusAL6PzB2Ywe9Igpg7vy7ThfckfmEm8WmVLjFJAdBeHNcJ7EWr2eY3wRl8MF9zvN8LLiXaVvUZrq2NHeS3rd1Wyfncl63dVUrS7kl0V9e3rjOifxtRhfbn1nBFMGZLFlCFZ9NF1maUbUUDEsoZqKP6zd/hIjfCiZn91Axv3VLFxTxUf761iw54qNu2toqaxBYD4OGNMTjpnjcpm0uA+TM7LYsqQPvRN09XVpHtTQMSa9kZ4hd65Cu2N8K7zBplHXaRGeCFobXXsrqynuLS6/ba5tJrifdWU1TS2r5ednsT43ExuLBjG+EGZTM7rw7jcTJ2QJj2SAiIWHNzhHTbasFCN8EJ2sLaRrftr2m9b9tewdZ/3uK6ppX29vmmJjM3JYPakXPJzM5kwKJNxuZkMyEjSiWjSa+i3TrS0NcIrKoTdq7xlORO98YQJc2HwGZp5dBJaWh17KuvZUVbL9gO1fFJWwycHatleVsu2/TVU1je3rxsfZwzrl8rIAemcPTqbMTkZjB3o3fqnKwhEFBCR4hzsXOENMhctVCO8k9TY3Mreynp2Hqxj18E6dpbXUVJex47yWkrKvWXNra59/fg4Y2i/VIZnp3Ht1DxGZKczakA6o3LSGdYvjaQEnRMicjQKiDB1bIRXtBCqdnmN8EZdoEZ4Adp++e+trGd3RT17KurZU+nd76rwfvmXVjXg3OGvy8lMZmi/VKYO68vc0wcztF8aw7JTGZGdTl7fFBLiFQIiJ0MB0dWa6ryrrBUt7NAIL9XrdTThu72uEZ5zjprGFvZVNbC/uoF9Vd6ttKqevZUNlFY1UFpZT2lVw2GDwW3SkuIZlJXC4KwULhqXQ17fVPKyUr37vink9U3VALFISBQQXaG9Ed4CvxFerd8I7ypvkHnMLEjqORdoqW9qoaymsf1WXtvI/upGDlQ3cKC6kQM1Dd7zGi8M2noKdRQfZ+RkJJPbJ5mh/dKYPqIfAzNTGJSVzKCsVAZnpTAoK4XM5ASNBYhEiQLiZFXtOTTzaOtiaG32G+Hd3G0a4dU3tVBR18TB2iYq6g7dDtY2crC2iYN1jZTXNlFR20S5v6yspvGw2T4dJcQZ2elJ9M9IZkBGEiP6p5GTkUxOpncbkHHovn96EnE6g1gkpikgTkRbI7yiQihZSnsjvHPu9gaZh0yPSCM85xwNza1U1TdT3dBMTUMzVfXNVNU3td9X+8sq/eeV9c1U1jVRWd9EZV0zlfVNNDZ/+i/7NnEGfdOS6JuaSN+0RHL7pDB+UCbZaUlkZySRnZZEv/Qk+qd79wPSk+mTqr/2RXoSBcSxOAd71hwaZD7BRnhtv8jrm1qobfRuhx43U+cvq21qobahuX1527o1Dc3UNDZT3eB93XveQnVDMy2tLvAzO0pKiKNPSgJ9UhLJTE2kT0oCQ/ql0ifFe9wnNZEsPwCy2h6nJpGVmkhmSoL+whfp5Xp9QBysbWTjnioamltpaG6lsamJjNLl5O56nSF7FpFZtxOHUZI5lfVD72V15vnssVzqd7bQsK2G+qYPqW9qob65hfqmVur8EKjzb0fOuDmelMQ40pMSSE2KJyM5gfTkBLJSE8nLSiE9OYH0pHgyUhLISE707+PJSE4kPTneC4KUBDKSE8hISdCFZETklMRcQJjZlcBPgXjgv5xzD4f5eSu2l/PFJ97n3Lh1zI5byuXxy8mxShpcAu+2TuHV1jm83nImFY19SamMJznBSE7YT0pSPCkJ8aQkxpGaFE9WaiIpifGkJMaTmhRHamI8qYnxpCT594nxpPmP05IS/HUSSEuKJy3ZX5YYr86eIhIzYiogzCwe+DlwOVACLDWzBc659V3+YX4jvPPWLmB95mskNFXTkphOzYhZ7BtzFa1jLmdqehYzE+N4KD5Oc+lFpNeJqYAAZgDFzrktAGb2LHAd0LUBse4P8Pu7oKWB5LT+MOV6mHgN8aMuoo8a4YmIALEXEEOAHR2elwBnd1zBzOYD8/2n1Wa28SQ/awCwHyrxdlp+fpJv063429yraJt7B23ziRnRmZViLSCCDsAfNszrnHsEeOSUP8hsmXOu4FTfpzvRNvcO2ubeIRLbHGsH1kuAYR2eDwV2RakWEZFeLdYCYimQb2ajzCwJuAlYEOWaRER6pZg6xOScazaze4BX8aa5Pu7GsrBqAAAGG0lEQVScWxfSx53yYapuSNvcO2ibe4fQt9nciZ7JJSIivUKsHWISEZEYoYAQEZFAPT4gzOxKM9toZsVm9kDA15PN7Dn/6x+Y2cjIV9m1OrHNXzOz9Wa22swWmVmn5kTHsuNtc4f1bjAzZ2bdfkpkZ7bZzD7rf6/Xmdn/RLrGrtaJn+3hZvamma30f77nRKPOrmJmj5tZqZmtPcrXzcx+5v97rDazM7u0AOdcj73hDXRvBkYDScBHwKQj1vk74Jf+45uA56JddwS2+RIgzX/8pd6wzf56mcBiYAlQEO26I/B9zgdWAv385wOjXXcEtvkR4Ev+40nAtmjXfYrbfCFwJrD2KF+fA7yMdw7ZTOCDrvz8nr4H0d66wznXCLS17ujoOuBJ//FvgVnWvS9qcNxtds696Zyr9Z8uwTvfpDvrzPcZ4AfA/wPqI1lcSDqzzV8Afu6cKwdwzpVGuMau1pltdkAf/3EW3fw8KufcYqDsGKtcBzzlPEuAvmY2uKs+v6cHRFDrjiFHW8c51wxUAP0jUl04OrPNHc3D+wukOzvuNpvZNGCYc25hJAsLUWe+z+OAcWb2rpkt8Tsld2ed2ebvA39rZiXAS8CXI1Na1Jzo//cTElPnQYTguK07OrlOd9Lp7TGzvwUKgItCrSh8x9xmM4sD/g24PVIFRUBnvs8JeIeZLsbbS3zbzKY45w6GXFtYOrPNNwNPOOd+bGbnAL/xt/nol0/s3kL9/dXT9yA607qjfR0zS8DbLT3WLl2s61S7EjO7DPgWcK1zriFCtYXleNucCUwB3jKzbXjHahd084Hqzv5s/8k51+Sc2wpsxAuM7qoz2zwPeB7AOfc+kILX1K6nCrU9UU8PiM607lgA3OY/vgF4w/mjP93UcbfZP9zyK7xw6O7HpeE42+ycq3DODXDOjXTOjcQbd7nWObcsOuV2ic78bP8Rb0ICZjYA75DTlohW2bU6s83bgVkAZjYRLyD2RbTKyFoA3OrPZpoJVDjndnfVm/foQ0zuKK07zOwfgWXOuQXAY3i7ocV4ew43Ra/iU9fJbf4RkAG84I/Hb3fOXRu1ok9RJ7e5R+nkNr8KzDaz9UAL8HXn3IHoVX1qOrnN9wOPmtlX8Q613N6d/+Azs2fwDhEO8MdVvgckAjjnfok3zjIHKAZqgTu69PO78b+diIiEqKcfYhIRkZOkgBARkUAKCBERCaSAEBGRQAoIEREJpICQXsHM3jKzK45Ydp+Z/edxXlfdRZ9f7d/nmdlvOyx/xu/C+VUzu93M8o7xHv9uZhf6j99qO9HPzF43s35dUadIRwoI6S2e4dPnuNzkL48Y59wu59wNAGY2CDjXOXe6c66tFUhgQJhZNjDTb952pN/gdSUW6VIKCOktfgvMNbNkAP+6H3nAO2aW4V8XY4WZrTGzT3WCNe+aIXM6PH/CzP7azOLN7EdmttTfE7jrWEWY2cgOvf1fAwaa2Soz+w5eX6yn/eepR7z0BuCVo7ztArweRCJdSgEhvYJ/BvGHQFtH07Zrfzi89t/XO+fOxGtN8eOAlu/PAp8D8Ns8zMI7i3UeXnuDs4CzgC+Y2ahOlnUtsNk5N9U59wNgGfA3/vO6I9Y9D1h+lG0rB5LNrDt3IZYYpICQ3qTjYaaOh5cM+KGZrQZex2uXnHvEa18GLvX3QK4CFvu/xGfj9cJZBXyA1yo+jIZ4gzl2T6FSjnJ4SuRk9eheTCJH+CPwE/+yjKnOuRX+8r8BcoDpzrkmv+NrSscXOufqzewt4Aq8PYmO4fJl59yrIdded2RNR0jx1xHpMtqDkF7DOVcNvAU8zuGD01lAqR8OlwBHu0b3s3jN0C7AaxiHf/8lM0sEMLNxZpZ+kiVW4bUmD1IEjA36gn84bBCw7SQ/VySQAkJ6m2eAM/B+2bd5Gigws2V4exMbjvLa1/CuEfy6f8lLgP8C1gMr/MHnX3Hye+ZPAL88yiD1i3hdPYNMB5b4V0QU6TLq5irSTZjZO8DcI68IZ2Y/BRY45xZFpzLpqbQHIdJ93A8MD1i+VuEgYdAehIiIBNIehIiIBFJAiIhIIAWEiIgEUkCIiEggBYSIiAT6X7qYwUyaMJt2AAAAAElFTkSuQmCC\n", "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "curve(50, 20, 'eqperc')" ] }, { "cell_type": "code", "execution_count": 8, "metadata": { "collapsed": false }, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYgAAAEKCAYAAAAIO8L1AAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMS4yLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvNQv5yAAAIABJREFUeJzt3Xd4VGXax/HvnYRUEgIJIaGEGpp0IiIIoqAgKKyKrlhRV3wtuxYsoLi6ulZU1rKrYlk7imKlFwVEpEsPHQKBUNJ7nef94wxuxAECZObMTO7PdXFlZnImc5+U+XHO85z7EWMMSiml1LEC7C5AKaWUd9KAUEop5ZIGhFJKKZc0IJRSSrmkAaGUUsolDQillFIuuS0gRKSZiPwoIikisklE7nE+/oSI7BeRtc5/Q6s8Z7yI7BCRrSIy2F21KaWUOjlx13UQIpIAJBhj1ohIJLAa+BNwNVBgjHnxmO07AlOAXkBjYD7Q1hhT6ZYClVJKnZDbjiCMMenGmDXO2/lACtDkBE8ZAXxmjCk1xuwGdmCFhVJKKRsEeeJFRKQF0B1YDvQF7haRG4FVwFhjTDZWeCyr8rQ0XASKiIwBxgBERET0bN++vVtrV0opf7N69eoMY0zDk23n9oAQkbrANOBeY0yeiLwBPAUY58eXgFsAcfH0P5z/MsZMBiYDJCcnm1WrVrmrdKWU8ksiklqd7dw6i0lE6mCFwyfGmK8AjDGHjDGVxhgH8Db/O42UBjSr8vSmwAF31qeUUur43DmLSYB3gRRjzMtVHk+ostnlwEbn7e+Aa0QkRERaAknACnfVp5RS6sTceYqpL3ADsEFE1jofewQYJSLdsE4f7QFuBzDGbBKRqcBmoAK4S2cwKaWUfdwWEMaYJbgeV5h5guc8DTztrpqUUkpVn15JrZRSyiUNCKWUUi5pQCillHJJA0IppZRLGhBKKaVc0oBQSinlkgaEUkoplzQglFJKuaQBoZRSyiUNCKWUUi5pQCillHJJA0IppZRLGhBKKaVc0oBQSinlkgaEUkoplzQglFJKuaQBoZRSyiUNCKWUUi5pQCillHJJA0IppZRLGhBKKaVc0oBQSikfUekwZBeWsW/bWo4cOuD21wty+ysopZT6g7IKB9lFZWQWlJFdVEZWYRk5RWVkF5X/7nZOURk5xeVkF5ZRWlLEnUHf8n+B37MpfgQN73jPrTVqQCilVA2oqHSQVVRGRn4ZmYWlZBaUkVlYRmbB0dulZBwNg4Iy8ksrjvu1IkOCiI6oQ/3wYKLDg2kRG0GP8rUM3/8S9Uv2sa/pMKIunOD2fdKAUEqp4zDGkFdcweH8Eo7kl3I4v9T50bp/pKCUjPwyMgpKySoqw5g/fo2gAKFBRDAxdUOIiQgmsUE4DSKCf/sXExFMfeft6PA6RIcFExxU5ex/wWGY8yhsmQoNWsFVX9Os9YUe2X8NCKVUrVRSXsnhvFLSc4s5mFfCwdwSDuWVcii/hMN5ztt5JZRWOP7w3JCgAOKiQoitG0JiTDg9mtenYWQIDesGE1s3xAqDusHERoQQFRaEiJx6gQ4HrPkA5j8OZUXQ/yHodz/UCauBva8eDQillN8pr3RwMLeEAznFpOeWcCC3mPSc/90/mFdCVmHZH54XHhxIfFQocVEhdE+MplFUKHGRIcRFhdKwbggNI0OIiwohMuQ03/Sr69AmmH4f7FsOzfvCpZOgYTv3vd5xaEAopXxOWYWD/TnF7MsqIi27mP05zo/ZxezPKeZQXgmOY073RIfXIaFeGI3rhdI9MZqEeqE0igoloV4Y8fVCaBQVSmRoHXt26KiyIlj0PPzyOoREwYj/QLdrwZ1hdAIaEEopr2OMIbuonNTMQvZmFbEno4i9WUXsyy4iLauI9LyS353vDwoQ4uuF0iQ6jHNbx9A0Oowm9cNoHB1mhUJ0KOHBXv52t20uzBwLOXuh2/Vw0ZMQEWNrSV7+HVNK+bPconJ2ZRSwJ7OQ3UcK2ZVRyJ7MQlIzi8gv+f0sn0ZRISQ2CKd3qxiaNggnsUE4zeqH0bRBOPFRoQQG2PO/7DOWlw6zH4bN30JsWxg9A1qcZ3dVgAaEUsrNKh2GtOwidhwuYOeRAnYeLmTnkQJ2ZRT+bhwgQKBZg3Cax0TQI7E+ic7bLWLCadYgnNA6gTbuhRs4KmHlO7DgKagsgwsmQN+/QVCI3ZX9RgNCKVUjKh2G1MxCth0qYPuhfLYdLmDbwXx2ZxRSVvm/mUCxdUNo1TCCwWfF0yo2gpaxEbSIjSCxQfjvp3f6swNrYfq9cOBXaHUBDHsJYlrbXdUfaEAopU5ZRkEpW9Lz2XIwj83peWw9mM+OwwW/mxLatH4YbRtFMqBdQ1o3rEvruLq0bhhBdHiwjZXbrDQffnwGlr8J4TFw5bvQ6UrbBqFPRgNCKXVcDochNauIjftz2XQgj00HcklJzyejoPS3beIiQ2ifEEWf1jEkNYqkXaNI2sTVJSJE315+YwxsmQ6zHoa8/dDzZhj0OITVt7uyE9KfoFIKsMJgd2Yh69NyWJ+Wy6b91tFBgbMlRJ1A+e2IoH18JB0TomgXH0lMXe85Z+6VcvbCzIdg2yyIOwuueh+a9bK7qmpxW0CISDPgQyAecACTjTGviEgD4HOgBbAHuNoYky3WVSevAEOBImC0MWaNu+pTqrY7mFvC2n3ZrEvLZd2+HDbsz/1t5lBYnUA6JERyRY8mdGpcj7OaRJEUF1l7xghqQmU5LHsDFj5r3b/oSeh9JwTafK3FKXDnEUQFMNYYs0ZEIoHVIjIPGA0sMMY8JyLjgHHAw8AlQJLz3znAG86PSqkzVFbhYHN6HmtSs1mzN5s1qdkcyC0BrCOD9vFRDO/amK7NounaNJo2cXV9d9qoN9i30hqEPrQRkgbDsBchOtHuqk6Z2wLCGJMOpDtv54tICtAEGAEMcG72AbAQKyBGAB8aYwywTESiRSTB+XWUUqegoLSCNanZrNidxYrdWaxLy/ltALlxvVB6NK/PXxLr0z0xmg4JUf43hdQuxTmw4ElY9R5ExsPVH0KH4V47CH0yHhmDEJEWQHdgOdDo6Ju+MSZdROKcmzUB9lV5Wprzsd8FhIiMAcYAJCb6XiIr5Q65xeXOMMhkxe4sNh7Io9JhCAwQOjWO4rpzmtOzeX16NI8moZ7nmr3VGsbAxmkw5xEoPALn3A4XPAqhUXZXdkbcHhAiUheYBtxrjMk7QYMrV5/4Q/NcY8xkYDJAcnKyi+a6Svm/4rJKVqdm8/PODJbuzGRDWg4OA8FBAXRrFs2dA1rTq2UDuifWp67OJnKvrF0wYyzs/AESusG1n0Pj7nZXVSPc+psjInWwwuETY8xXzocPHT11JCIJwGHn42lAsypPbwq4f009pXyAMYbN6Xks3pbBom2HWZOaQ1mlg6AAoVuzaO6+MIk+rWPonhhNSJCeLvKIijJY+iosnggBdWDI89DrNgjwn++/O2cxCfAukGKMebnKp74DbgKec378tsrjd4vIZ1iD07k6/qBqs+zCMn7akcGirUdYvP0IR/Ktaw/ax0dyU5/m9GkTS68WDfR6AzukLrXacR/ZYo0xXPI8RDW2u6oa587frL7ADcAGEVnrfOwRrGCYKiK3AnuBq5yfm4k1xXUH1jTXm91Ym1JexxjDziMFzE85zPzNh1izNxuHsdpU90tqyPltG9I/KZa4qFC7S629irJg3mPw68dQrxmM+hzaDbG7Krdx5yymJbgeVwAY6GJ7A9zlrnqU8kYVlQ5W7slmfsohFqQcYk9mEQCdmkTx1wuTGNCuIV2aRuuUU7sZA+umwNwJ1kylPn+DAeMgOMLuytxKj02V8rDySgdLd2Yya0M6czcfIquwjOCgAPq0juEv/VoxsEOczjTyJke2wYz7Yc9P0PRsuPRfEN/J7qo8QgNCKQ8orajk5x0ZzNxwkHmbD5FbXE7dkCAGdohjyFnx9G/bUMcSvE15CSx5GZZMstaBvnQS9BgNAbXnanL9jVTKTRwOw4o9WXy7dj8zNxwkt7icyNAgLurYiKGdEjgvKVYvUPNWuxbC9Pshayd0vgoGPwN14076NH+jAaFUDUtJz+Obtfv5fu0BDuSWEB4cyMUdGzG8W2POa9NQ+xl5s4Ij1sVuG6ZC/ZZw/VfQ5g9DprWGBoRSNSC7sIxv1u5n6qo0UtLzCAwQ+ifF8vAl7bmoYyPvXw+5tnM4YM0HMP9xKCuC/g9Cv7HWqaVaTH9rlTpNlQ7Dkh0ZTF21j3mbDlFW6aBL03o8OeIshnVO0DbYvuLQJuuahn3LoXlfa6yhYTu7q/IKGhBKnaKDuSV8umIvX67ax4HcEqLD63Bd70SuTm5GhwTf7r1Tq5QVwaLn4ZfXISQKRvwHul3rs4313EEDQqlqMMbwy65MPvollbmbD+Ewhn5JDXl0WEcGdYzT9ha+ZttcmDnWWsyn23Vw0VMQEWN3VV5HA0KpEygoreCrNWl89Esq2w8XEB1eh1vPa8l15yTSPMa/L5LyS3npMPth2PwtxLaF0TOgxXl2V+W1NCCUcuFATjH//Xk3n63YR35pBV2a1mPiyC5c1rWxTk31RY5KWPmutVZDZRlcMAH6/g2CdJzoRDQglKpi4/5c3vlpF9PXp2OAYZ0TuLlvC7onevfi8uoE0tfB9/fAgV+h1QUw7CWIaW13VT5BA0LVesYYFm/PYPLinfy8I5OI4EBu6tOCm/u2oGn9cLvLU6erNB9+fAaWvwnhMXDlu9DpSh2EPgUaEKrWMsYwP+Uwry7Yzob9uTSKCmHcJe0Z1SuRemG+s7C8ciFlOsx6CPL2Q8+bYdDjEKZHgadKA0LVOg6HYe7mQ7y6YDub0/NIbBDO81d25vLuTfUqZ1+Xs88Khq0zIe4suOp9aNbL7qp8lgaEqjUcDsPsTQd5dcF2thzMp2VsBC9e1ZUR3RpTJ1CDwadVVsDyN+DHZwEDFz0Jve+EQD0SPBMaEMrvHR1jeGH2FjYdyKN1wwj+9eduXNolgSANBt+Xtgq+vxcObYCkwTB0ItRvbndVfkEDQvm1dftyeH72FpbuzKRZgzAm/bkrw7s20QV4/EFJrjVtdeW7EBkPV39oLf+pg9A1RgNC+aVdRwp4ce5WZm44SExEME9c1pFrz2muYwz+wBjY9BXMHg+FR+Cc2+GCRyFU25zUNA0I5VdyisqYNG8bHy/fS0hQAPcMTOK2/q2oq4vx+Ies3TBjLOxcAAldYdRn0KSH3VX5Lf2rUX6h0mH4dMVeXp67ldzicq49J5F7BralYaReKesXKspg6auweCIEBMGQ5+Ds2yBQ38LcSb+7yuct35XJE99vJiU9j96tGvD4ZWdpV1V/krrUasd9ZAt0uAyGPA/1mthdVa2gAaF8VnpuMU/PSGH6+nSaRIfxn+t6cEmneEQHKf1DURbM+zv8+hHUawajPod2Q+yuqlbRgFA+p9Jh+OiXPUycs5UKh+HeQUnc3r81YcHaRM8vGAPrPoO5j0JxDvT5GwwYB8HaPdfTNCCUT9lyMI9x0zawdl8O/ds25Ok/daJZA+2X5Dcytlunk/b8BE3PtlZ3i+9sd1W1lgaE8gkl5ZW8umA7kxfvol5YHV65phvDuzbW00n+orwElrwMSyZBUBgMe9nqoRSg05LtpAGhvN6yXZmMm7aePZlFjOzZlEeHdqB+RLDdZamasmshTL8fsnZCp5Ew+BmIbGR3VQoNCOXFSsormThnK+8u2U3zmHA++cs59G0Ta3dZqqYUHLHGGdZ/DvVbwvVfQZuBdlelqtCAUF5p4/5c7vt8LdsPF3BD7+aMH9qe8GD9dfULDgf8+iHMexzKCqH/g9BvLNQJs7sydQz9i1NepaLSwZuLdvKv+dtpEBHMB7f04vy2De0uS9WUQ5utQeh9y6B5X2sQumE7u6tSx6EBobzGnoxC7p+6ljV7cxjWJYGn/9SJ6HAda/ALZUWw+AVY+hqERMGIf0O367SxnpfTgFBe4du1+3nkqw0EBojOUPI32+dZ/ZNyUq1QuOgpiIixuypVDRoQylYl5ZU8OX0zny7fS8/m9XltVHcaR+u5aL+Qlw6zx8HmbyC2LYyeAS3Os7sqdQo0IJRtdmcUcucna0hJz+P281vxwMXtdGU3f+CotNZo+OEpqCiFCyZA379BkDZO9DUaEMoW3687wPivNhAUKLw3OpkL2+u8d7+Qvs5a3e3AGmg1wLrgLaa13VWp0+S2/66JyHsiclhENlZ57AkR2S8ia53/hlb53HgR2SEiW0VksLvqUvYqq3Dw2Dcb+euUX2nbqC4z/tZPw8EflBbA7Edg8gDI3QdXvAM3fKPh4OPceQTxPvA68OExj08yxrxY9QER6QhcA5wFNAbmi0hbY0ylG+tTHpZRUModH69m5Z5sbuvXkoeGtNdTSv4gZTrMegjy9lvtMQY9DmH17a5K1QC3BYQxZrGItKjm5iOAz4wxpcBuEdkB9AJ+cVN5ysM2HchlzIerySgo5dVR3RnetbHdJakzlbPPCoatMyHuLLjqfWjWy+6qVA2yYwzibhG5EVgFjDXGZANNgGVVtklzPvYHIjIGGAOQmJjo5lJVTZi5IZ2xU9cRHV6HL/+vD52b1rO7JHUmKitg+Rvw47NgHHDRk9D7TgisY3dlqoZ5+vj+DaA10A1IB15yPu5qwrtx9QWMMZONMcnGmOSGDfUKW2/mcBhenreNOz9ZQ4eESL69u6+Gg69LW2WNM8ydYE1ZvWs59L1Hw8FPefQIwhhz6OhtEXkbmO68mwY0q7JpU+CAB0tTNayorIL7P1/H7E0HGdmzKU9f3omQIF3Qx2eV5MKCJ63pq5HxcPWH0GG4Xgnt5zwaECKSYIxJd969HDg6w+k74FMReRlrkDoJWOHJ2lTNySwo5ZYPVrEhLYcJwzpw63kt9apoX2UMbPoKZo+HwiNwzu1wwaMQqmt+1wZuCwgRmQIMAGJFJA14HBggIt2wTh/tAW4HMMZsEpGpwGagArhLZzD5ptTMQm56bwXpuSW8eX1PLj4r3u6S1OnK2g0zH4Ad8yGhK4z6DJr0sLsq5UFijMtT/T4hOTnZrFq1yu4ylNP6tBxueX8lFQ7DuzedTc/mOtXRJ1WUwdJXYfFECAiCCyfA2bdBoF5X6y9EZLUxJvlk2+lPXNWIhVsPc+cna35r0d26YV27S1KnI/UXmH4vHNkCHS6DIc9DPZcTClUtoAGhztgXq/Yx7qsNtGsUyfs3n01cVKjdJalTVZQF8/4Ov34E9ZrBqM+h3RC7q1I204BQZ+Ttxbt4emYK/ZJi+c91PYgM1emOPsUYWPeZtfRncQ70+SsMGA/BEXZXpryABoQ6ba8t2M5L87YxrEsCk67uRnCQts3wKRnbYcb9sHsxND3bWt0tvrPdVSkvogGhTpkxhhfnbuXfP+7kih5NmDiyK4EBOo3VZ5SXwJJJsORlCAqzOq72vBkCNODV72lAqFNijOGfM1J4d8luRvVK5Ok/dSJAw8F37FpkrQmdtRM6jYTBz0CkdtNVrmlAqGpzOAx//24jHy/by+g+LXj8so56AZyvKDhijTOs/xzqt4Drv4I2A+2uSnk5DQhVLZUOw/iv1jN1VRr/d35rHh7STsPBFzgc1sykeX+HskLo9wD0fwDq6LKu6uQ0INRJORyGh6et58vVadw7KIl7BiZpOPiCwynW6m77lkHzvtYgdMN2dlelfIgGhDohYwyPf7eJL1encd+gttwzKMnuktTJlBXB4hdg6WsQEgkj/g3drtPGeuqUaUCo4zLG8NysLXy0LJXbz2/F3wa2sbskdTLb58GMsZCTCl2vhYufgohYu6tSPkoDQh3Xqwt28NbiXdx4bnPGDWmvp5W8Wf5BmD0ONn0NMUlw03Ro2c/uqpSP04BQLr29eBeT5m9jZM+mPHHZWRoO3spRCaves9ZqqCi1WnH3vQeCQuyuTPkBDQj1Bx8tS+XpmSkM65LA81d20escvFX6OmsQ+sAaaDXAuuAtprXdVSk/ogGhfufrX9N47JuNDGwfx6Sru+kV0t6otAB+fMZaFzo8Bq54BzqP1EFoVeOqFRAicv+JPm+MeblmylF2WrztCA9+sZ4+rWP493U9tLeSN9oyA2Y+CHn7oedoGPQEhOm6G8o9qnsEkQycjbU0KMBlwGJgnzuKUp636UAud3y8mjZxdXnzhp6E1tH1o71KbhrMfAi2zoC4jjDyv5B4jt1VKT9X3YCIBXoYY/IBROQJ4AtjzF/cVZjynLTsIm7+70qiwurw/s29iNKW3d6jsgKWv2mdUjIOGPQPOPcuCNSfkXK/6gZEIlBW5X4Z0KLGq1Eel1tUzuj/rqS4vJIv/68P8fV0sR+vkbYapt8DBzdA0mAYOhHqN7e7KlWLVDcgPgJWiMjXgAEuBz5wW1XKI0rKK7ntw1XszSzig1t60S4+0u6SFEBJLix4Cla+A5HxcPWH0GG4DkIrj6tWQBhjnhaRWcDRK29uNsb86r6ylLs5HIaxU9exYk8Wr47qzrmtY+wuSRljXeg2ezwUHIJeY+DCCRAaZXdlqpaq9jRXY8waYI0ba1Ee9OysFGZsSOeRoe0Z3rWx3eWorN0w8wHYMR8SusKoKdCkh91VqVpOr4Oohaau2sfbP+3mxnObc1u/VnaXU7tVlMEvr8GiFyAgCIY8B2ffBoH6p6nsp7+Ftczq1GwmfL2R89rE8vdLdcEfW6X+Yq3udiQFOlwGQ56Hek3srkqp32hA1CLpucXc/tFqEqJDef3a7gQF6oVwtijKgvmPw5oPoV4zGPUZtLvE7qqU+gMNiFqiuKySMR+upqS8kim3nUN0eLDdJdU+xlhLfs55FIqzoc9f4fxxEFLX7sqUckkDohYwxvDQtPVsPJDL2zckk9RIp7N6XMYOmHEf7F4MTZLhxm8gvrPdVSl1QhoQtcB/Fu7k+3UHeGhIOwZ1bGR3ObVLeQksmQRLXoagMKvjas+bIUBP7ynvpwHh5+ZvPsSLc7cyvGtj7jhfW0F71K5FMON+yNwBnUbC4GcgUgNa+Q4NCD+2O6OQez9fS6fG9XhhZBedseQphRnWOMP6z6B+C7j+K2gz0O6qlDplGhB+qriskjs+Xk2dQNHurJ7icMCvH8G8v0NZIfR7APo/AHXC7K5MqdOiAeGn/v7tRrYeyue/o8+mSbS+Qbnd4RTrmoa9v0BiH7h0EsS1t7sqpc6IBoQfmrpyH1+sTuNvF7ZhQLs4u8vxb2VFsPgFWPoahETC8Neh23U6CK38ggaEn9l8II/Hvt1I3zYx3DOord3l+Lft861B6JxU6HotXPwURMTaXZVSNUYDwo/klZRz5yeriQ6vwyvXdNf1pN0l/yDMHmd1Xo1JgpumQ8t+J3+eUj7GbcfBIvKeiBwWkY1VHmsgIvNEZLvzY33n4yIir4rIDhFZLyLaxvIUGWN4+Mv17Msu5vVrexBbN8TukvyPoxJWvA2vnw1bZsIFj8IdP2s4KL/lzhOl7wNDjnlsHLDAGJMELHDeB7gESHL+GwO84ca6/NJ/f97DrI0HeXhIO85u0cDucvxP+np49yKrJXfj7nDnL3D+QxCkQaz8l9tOMRljFotIi2MeHgEMcN7+AFgIPOx8/ENjjAGWiUi0iCQYY9LdVZ8/Wbcvh2dmpnBRx0bavrumlRbAwmdh2RsQ3gCueBs6X6Wru6lawdNjEI2OvukbY9JF5OgUmybAvirbpTkf+0NAiMgYrKMMEhMT3VutDygsreDez9cSFxnCiyO76sVwNWnLDJj5EOSlQc/RMOgJCKtvc1FKeY63DFK7elczrjY0xkwGJgMkJye73KY2+cf3m9iTWciU23pTL7yO3eX4h9w0Kxi2zoC4jjByLiSeY3dVSnmcpwPi0NFTRyKSABx2Pp4GNKuyXVPggIdr8zkzN6QzdVUad13Qmt6tdE3pM1ZZASvegh+eBuOwjhjOvRsCNXhV7eTpq3m+A25y3r4J+LbK4zc6ZzP1BnJ1/OHEDuQUM27aero2i+Zevd7hzKWthrcHwJxHoEVfuGsZnHefhoOq1dx2BCEiU7AGpGNFJA14HHgOmCoitwJ7gaucm88EhgI7gCLgZnfV5Q8qHYb7Pl9LpcPwyp+7UUdXhjt9Jbmw4ClY+Q5ExsNVH0DHEToIrRTuncU06jif+kNbS+fspbvcVYu/eXPRTpbvzmLiyC60iI2wuxzfZAxs/gZmjYOCQ9BrDFw4AUKj7K5MKa/hLYPUqprW7cth0rxtDOucwMieTe0uxzdl74EZD8COeZDQFUZNgSZ6baZSx9KA8CGFpRXc89mvxEWG8MzlnXVK66mqLLea6i16AQICYfCz1pFDoP4ZKOWK/mX4kGdnpZCaVaRTWk/H3mXw/b1wJAXaXwqXvAD1mthdlVJeTQPCRyzZnsHHy/Zy63ktdUrrqSjKgvmPw5oPoV4zGPUZtLvE7qqU8gkaED4gv6Sch6etp1VsBA8Obmd3Ob7BGFg/1Zq2WpwNff4K54+DkLp2V6aUz9CA8AHPzEwhPbeYL+/oo0uHVkfGDmudht2LoEky3PgNxHe2uyqlfI4GhJdbtO0IU1bs4/bzW9EjUfsAnVBFKSyZBD+9BEFhMOwl6HmzNSCtlDplGhBeLLe4nIe/XE+buLrcp1dLn9juxdaa0Jk7oNOV1gylyEZ2V6WUT9OA8GL/nL6ZIwWlvHVDTz21dDyFGTDnUVj/GdRvAddPgzaD7K5KKb+gAeGlfthyiC9Wp3H3BW3o2iza7nK8j8MBaz+GuY9BWSH0ewD6PwB1wuyuTCm/oQHhhXKLyhk3bQPt4yP568A2dpfjfQ6nWKeT9v4CiX3g0kkQ197uqpTyOxoQXuifMzaTWVjGe6PPJiRITy39pqwIFk+Epa9CSCQMfx26XQcB2qxQKXfQgPAyS3dk8MXqNO4Y0JpOTerZXY732D7fmrqakwpdr4WLn4KIWLurUsqvaUB4kZLySsZ/vYEWMeHcMzDJ7nK8Q/5BmD0eNn0FMUlw03Ro2c/uqpSqFTQgvMgrC7aTmlnEp385R2ctOSph1Xuw4Enr+oYLHoW+90BQiN2VKVVraEB4iZT0PCYv3sXInk3p06aWnzpJXw/T74X9q6Hl+daAOu8xAAASA0lEQVQgdExru6tSqtbRgPAClQ7DuGnriQ6rw6NDO9hdjn1KC2Dhs7DsDQhvAFe8DZ2v0tXdlLKJBoQX+GDpHtal5fLKNd2oHxFsdzn22DIDZj4EeWnQczQMegLCtLWIUnbSgLBZWnYRL87dyoB2DRnetbHd5XhebhrMehi2TIe4jjByDiT2trsqpRQaELYyxvDYNxsxBp4a0al2rRBXWQEr3oIfn7EGpAc9AefeDYG6EJJS3kIDwkbT16fz49YjTBjWgWYNwu0ux3P2r7ZWdzu4HpIuhqETrT5KSimvogFhk/yScp6avplOTaK4uW9Lu8vxjJJcWPAUrHwH6jaCqz6AjiN0EFopL6UBYZNJ87ZzpKCUyTcmExjg52+QxsDmb2DWOCg4BL1ugwsnQKheKa6UN9OAsMHmA3m8v3Q3o3ol0s3fO7Vm74EZD8COeRDfBUZ9Ck162l2VUqoaNCA8zOEwPPbtRqLDg3nIn9eXriyHpa/BohesFd0GPwu9xkCg/sop5Sv0r9XDvlyTxurUbF4Y2YXocD+95mHvMmsQ+kgKtL8ULnkB6jWxuyql1CnSgPCgnKIynpu1hZ7N6zOyR1O7y6l5RVkw/wlY8wFENYVrpkD7oXZXpZQ6TRoQHjRxzlZyisp4asQ5BPjTwLQxsH4qzHkEirOt6xkGjIeQunZXppQ6AxoQHrJuXw6frtjLzX1a0rFxlN3l1JyMHdY6DbsXQZNkuPEbiO9sd1VKqRqgAeEBlQ7DhG820rBuCPdd5CfrPFSUwpJ/wU8vQVAoDHsJet5sDUgrpfyCBoQHfLpiLxv2W834IkP9oJXE7sUw/X7I3A6drrRmKEU2srsqpVQN04Bws+zCMl6cs5VzW8X4fjO+wgyYOwHWTbFaY1w/DdoMsrsqpZSbaEC42cvztlFQWsETw8/y3WZ8Dges/Rjm/d1as6HfWOj/INQJs7sypZQbaUC4UUp6Hp8sT+WG3s1pFx9pdzmn5/AWmH4f7F0KiX2s1d3i2ttdlVLKA2wJCBHZA+QDlUCFMSZZRBoAnwMtgD3A1caYbDvqqwnGGP7x/Saiwupw30Vt7S7n1JUXw+KJ8POr1nTV4a9Dt+sgIMDuypRSHmLnX/sFxphuxphk5/1xwAJjTBKwwHnfZ83aeJBlu7IYe1Fb37tiesd8+E9va4ZS56vg7lXQ4wYNB6VqGW86xTQCGOC8/QGwEHjYrmLOREl5JU/PSKF9fCSjeiXaXU715R+E2eNh01cQ0wZu+h5a9re7KqWUTewKCAPMFREDvGWMmQw0MsakAxhj0kUkztUTRWQMMAYgMdE733zfWrSL/TnFTLmtN0GBPvC/bocDVr8H8/9hXd8w4BE4714ICrG7MqWUjewKiL7GmAPOEJgnIluq+0RnmEwGSE5ONu4q8HTtzynmjUU7GNo5nnNbx9hdzskd3GA11tu/Clqebw1Cx7S2uyqllBewJSCMMQecHw+LyNdAL+CQiCQ4jx4SgMN21Hamnp2ZgjHwyNAOdpdyYqUFsPBZWPYGhDeAK962xht8dSquUqrGefz8h4hEiEjk0dvAxcBG4DvgJudmNwHferq2M7V8VybT16dz+/mtaVrfi9eY3jIT/n0O/PI6dL8e7loBXa7WcFBK/Y4dRxCNgK+dF40FAZ8aY2aLyEpgqojcCuwFrrKhttNW6TA8OX0zjeuFcsf5XnqKJjcNZj0MW6ZDXEcYOQcSe9tdlVLKS3k8IIwxu4CuLh7PBAZ6up6aMm1NGpsO5PHKNd0IC/ayhnWVFbBiMvz4NDgqYdATVkvuQD/oC6WUchtvmubqswpLK3hxzla6J0Z7X7+l/autQeiD6yHpYhg60eqjpJRSJ6EBUQPeWryLw/mlvHF9T+/pt1SSCz/8E1a8DXUbwVUfQMcROs6glKo2DYgzlJ5bzOTFO7msa2N6Nq9vdznW6m6bv4FZ46DgEPS6DS6cAKH17K5MKeVjNCDO0MTZW3EYeGhwO7tLgexUmPkAbJ8L8V1g1KfQpKfdVSmlfJQGxBlYn5bDV7/u544BrWnWwMZprZXl1pTVhc+DBFgL+PQaA4H641VKnT59BzlNxhj+OT2F2LrB3DnAxmmte5fD9Hvh8GZofylc8jzUa2pfPUopv6EBcZpmbzzIij1ZPHN5Z3uWES3KgvlPwJoPIKopXDMF2g/1fB1KKb+lAXEaSisqeXbWFto1iuTqZA//b90YWD8V5jwCxdnW9QwDxltrNiilVA3SgDgNHy5NZW9WER/e0suz3Vozd1qru+1eZA0+3/A1JHTx3OsrpWoVDYhTlF1Yxqs/bOeCdg3p37ahZ160ohSW/MtawCcoBIa+CMm3QICXXbGtlPIrGhCn6NUftlNYWsF4T3Vr3f2TddSQuR3OugKGPAuR8Z55baVUraYBcQpSMwv5eFkqfz67GW0bRbr3xQozYO4EWDfFao1x/TRoM8i9r6mUUlVoQJyCF2ZvpU5gAPcNauu+F3E4YO3HMO/vUJoP/cZC/wehTpj7XlMppVzQgKimNXuzmbEhnXsGJhEXFeqeFzm8xTqdtHcpJJ5rre4W5+ULDyml/JYGRDUYY3hmRgqxdUMY079Vzb9AeTEsngg/v2pNVx3+GnS7HgJ8YD1rpZTf0oCohrmbD7EqNZunL+9EREgNf8t2zIcZYyF7D3QdBRf/EyJia/Y1lFLqNGhAnER5pYPnZ22hTVxd/pzcrOa+cP4hmDMeNk6DmDZw0/fQsn/NfX2llDpDGhAn8dmKvezKKOSdG5Nr5qI4hwNWvwfzn4SKYhjwCJx3r3V9g1JKeRENiBPILynnX/O3c07LBgzsEHfmX/DgBmt1t/2rrKOFYZMgts2Zf12llHIDDYgTeGvRLjILy/jvsA5ntlJcaQEsfBaWvQFh9eHyydDlal3dTSnl1TQgjuNgbgnvLNnF8K6N6dI0+vS/0NZZMPNByN0HPW6CQU9AeIOaKlMppdxGA+I4Js3bRqXD8ODprhSXux9mPQRbpkNcR7hlDiT2rtkilVLKjTQgXNh+KJ8vVu9jdJ+Wp75SXGUFrJgMPz4NjkrriOHcuyHQhjUjlFLqDGhAuPD87C1EBAdx94WnOIC8f421ulv6OmhzEQx70eqjpJRSPkgD4hgrdmcxP+UwDw5uR4OI4Oo9qSQPfvgnrHwbIuLgqveh4590EFop5dM0IKowxvDsrBTio0K5pW/L6jwBNn8Ls8dB/kHodRtcOAFC67m/WKWUcjMNiCrmbDrIr3tzeP7KzoQFn2QxnuxUmPkAbJ8L8V3gmk+sVd6UUspPaEA4lVc6eGH2VpLi6nJljxOsM11ZDr/8GxY+BxIAg5+BXrdDoH4rlVL+Rd/VnD5fue/kLTX2LrfacR/eBO0vhUueh3onCBOllPJhGhBAYWkF/5q/nV4tjtNSozgb5j8Bq9+HqKZwzRRoP9TTZSqllEdpQADv/LSbjIJSJt/Y8/ctNYyBDV/AnEegKMu6nmHAeGvNBqWU8nO1PiAyCkqZvHgnl3SKp0di/f99InMnzLgfdi20Bp+v/woSuthWp1JKeVqtD4gN+3MJDJD/tdSoKIWfX4HFL1otuIe+CMm3QMBJZjUppZSfqfUBcUG7OJY/Msia1rr7J2sQOnM7nHUFDHkWIuPtLlEppWxR6wMCIKw8B2ZMgHWfQnRzuG4aJA2yuyyllLKV1wWEiAwBXgECgXeMMc+57cWMgV8/hnmPQWk+nHc/9H8Qgk+xQZ9SSvkhrwoIEQkE/g1cBKQBK0XkO2PM5hp/sSNbrdNJqT9D4rlw6SSI61DjL6OUUr7KqwIC6AXsMMbsAhCRz4ARQM0GxKZvYNpfIDgChr8G3a6HgBpYb1oppfyItwVEE2BflftpwDlVNxCRMcAY590CEdl6mq8VC5kZjL8JuOk0v4TPiQUy7C7Cw3Sfawfd51PTvDobeVtAuOqPbX53x5jJwOQzfiGRVcaY5DP9Or5E97l20H2uHTyxz952XiUNaFblflPggE21KKVUreZtAbESSBKRliISDFwDfGdzTUopVSt51SkmY0yFiNwNzMGa5vqeMWaTm17ujE9T+SDd59pB97l2cPs+izHm5FsppZSqdbztFJNSSikvoQGhlFLKJb8PCBEZIiJbRWSHiIxz8fkQEfnc+fnlItLC81XWrGrs8/0isllE1ovIAhGp1pxob3ayfa6y3UgRMSLi81Miq7PPInK182e9SUQ+9XSNNa0av9uJIvKjiPzq/P326ZW9ROQ9ETksIhuP83kRkVed34/1ItKjRgswxvjtP6yB7p1AKyAYWAd0PGabO4E3nbevAT63u24P7PMFQLjz9h21YZ+d20UCi4FlQLLddXvg55wE/ArUd96Ps7tuD+zzZOAO5+2OwB676z7Dfe4P9AA2HufzQ4FZWNeQ9QaW1+Tr+/sRxG+tO4wxZcDR1h1VjQA+cN7+Ehgov1tWzuecdJ+NMT8aY4qcd5dhXW/iy6rzcwZ4CngBKPFkcW5SnX2+Dfi3MSYbwBhz2MM11rTq7LMBopy36+Hj11EZYxYDWSfYZATwobEsA6JFJKGmXt/fA8JV644mx9vGGFMB5AIxHqnOPaqzz1XdivU/EF920n0Wke5AM2PMdE8W5kbV+Tm3BdqKyM8isszZKdmXVWefnwCuF5E0YCbwV8+UZptT/Xs/JV51HYQbnLR1RzW38SXV3h8RuR5IBs53a0Xud8J9FpEAYBIw2lMFeUB1fs5BWKeZBmAdJf4kIp2MMTlurs1dqrPPo4D3jTEvici5wEfOfXa4vzxbuPX9y9+PIKrTuuO3bUQkCOuw9ESHdN6uWu1KRGQQ8Cgw3BhT6qHa3OVk+xwJdAIWisgerHO13/n4QHV1f7e/NcaUG2N2A1uxAsNXVWefbwWmAhhjfgFCsZra+Su3tify94CoTuuO7/hfO9eRwA/GOfrjo066z87TLW9hhYOvn5eGk+yzMSbXGBNrjGlhjGmBNe4y3Bizyp5ya0R1fre/wZqQgIjEYp1y2uXRKmtWdfZ5LzAQQEQ6YAXEEY9W6VnfATc6ZzP1BnKNMek19cX9+hSTOU7rDhF5ElhljPkOeBfrMHQH1pHDNfZVfOaquc8TgbrAF87x+L3GmOG2FX2GqrnPfqWa+zwHuFhENgOVwIPGmEz7qj4z1dznscDbInIf1qmW0b78Hz4RmYJ1ijDWOa7yOFAHwBjzJtY4y1BgB1AE3Fyjr+/D3zullFJu5O+nmJRSSp0mDQillFIuaUAopZRySQNCKaWUSxoQSimlXNKAULWCiCwUkcHHPHaviPznJM8rqKHXL3B+bCwiX1Z5fIqzC+d9IjJaRBqf4Gv8S0T6O28vPHqhn4jMF5H6NVGnUlVpQKjaYgp/vMblGufjHmOMOWCMGQkgIvFAH2NMF2PM0VYgLgNCRBoAvZ3N2471EVZXYqVqlAaEqi2+BC4VkRAA57ofjYElIlLXuS7GGhHZICJ/6AQr1pohQ6vcf19ErhSRQBGZKCIrnUcCt5+oCBFpUaW3/1wgTkTWishjWH2xPnHeDzvmqSOB2cf5st9h9SBSqkZpQKhawXkF8QrgaEfTo2t/GKz235cbY3pgtaZ4yUXL98+APwM42zwMxLqK9Vas9gZnA2cDt4lIy2qWNRzYaYzpZox5ClgFXOe8X3zMtn2B1cfZt2wgRER8uQux8kIaEKo2qXqaqerpJQGeEZH1wHysdsmNjnnuLOBC5xHIJcBi55v4xVi9cNYCy7FaxbujIV4CJ+4pdJjjnJ5S6nT5dS8mpY7xDfCyc1nGMGPMGufj1wENgZ7GmHJnx9fQqk80xpSIyEJgMNaRRNVw+asxZo6bay8+tqZjhDq3UarG6BGEqjWMMQXAQuA9fj84XQ847AyHC4DjrdH9GVYztH5YDeNwfrxDROoAiEhbEYk4zRLzsVqTu5ICtHH1CefpsHhgz2m+rlIuaUCo2mYK0BXrzf6oT4BkEVmFdTSx5TjPnYu1RvB855KXAO8Am4E1zsHntzj9I/P3gTePM0g9A6urpys9gWXOFRGVqjHazVUpHyEiS4BLj10RTkReAb4zxiywpzLlr/QIQinfMRZIdPH4Rg0H5Q56BKGUUsolPYJQSinlkgaEUkoplzQglFJKuaQBoZRSyiUNCKWUUi79PxL9YTsep1dzAAAAAElFTkSuQmCC\n", "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" }, { "data": { "text/plain": [ "" ] }, "execution_count": 8, "metadata": {}, "output_type": "execute_result" } ], "source": [ "interact(curve, \n", " R=(5., 100.), \n", " Ccv=(5., 200.), \n", " characteristic=['linear', 'eqperc'])" ] } ], "metadata": { "anaconda-cloud": {}, "kernelspec": { "display_name": "Python 3", "language": "python", "name": "python3" }, "language_info": { "codemirror_mode": { "name": "ipython", "version": 3 }, "file_extension": ".py", "mimetype": "text/x-python", "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython3", "version": "3.6.2" } }, "nbformat": 4, "nbformat_minor": 2 }